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1. INTRODUCTION 

Our goal in this paper is to derive, in terms of a triplet of constant matrices, explicit 
formulas for exact solutions to the sine-Gordon equation 

Uxt = smu, (1.1) 

where u is real valued and the subscripts denote the partial derivatives with respect to the 

spatial coordinate x and the temporal coordinate t. Under the transformation 

t t 

X ^ ax -\ — , t ^ ax , 

a a 

where a is a positive constant, (1.1) is transformed into the alternate form 

Uxx-utt = smu, (1.2) 

and hence our explicit formulas can easily be modified to obtain explicit solutions to (1.2) 
as well. Let us note that one can omit a multiple of 2tt from any solution to (1.1). We 
are interested in solutions to (1.1) where Ux{x,t) — > as a; — )■ ±oo for each fixed t, and 
hence without any loss of generality we will normalize our solutions so that u{x, t) — > as 

X +00. 

The sine-Gordon equation arises in applications as diverse as the description of surfaces 
of constant mean curvature [10,16], one-dimensional crystal dislocation theory [17,23,40,41], 
magnetic flux propagation in Josephson junctions (gaps between two superconductors) 
[29,31], condensation of charge density waves [11,22,36], wave propagation in ferromag- 
netic materials [19,27,30], excitation of phonon modes [35], and propagation of deforma- 
tions along the DNA double helix [18,26,38,43]. 

The literature on exact solutions to (1.1) is large, and we will mention only a few and 
refer the reader to those references and further ones therein. For a positive constant a, by 
substituting 

1 /U(ax + a~^t)\ 
«(.,0=4ta.-( ^|^^_^_,^j j, (1.3) 
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into (1.1) and solving the resulting partial differential equations for U and V, Steuerwald 
[42] has catalogued many exact solutions to the sine-Gordon equation in terms of elliptic 
functions. Some of these solutions, including the one-soliton solution, two-soliton solutions 
modeling a soliton-soliton and soliton-antisoliton collision, and the breather solution, can 
be written in terms of elementary functions [25,37], while the n-soliton solutions can be 
expressed as in (1.3) where U and V are certain determinants [34,39,45]. The same sepa- 
ration technique can also be used to find exact solutions to the sine-Gordon equation on 
finite (a;-|-t)-intervals [13]. Solutions to the sine-Gordon equation with initial data specified 
on invariant algebraic manifolds of conserved quantities can be written explicitly in terms 
of Jacobi theta functions [24]. The ordered exponential integrals appearing in such solu- 
tions can be evaluated explicitly [9,28]. Let us also mention that some exact solutions to 
the sine-Gordon equations can be obtained via the Darboux or Backlund transformations 
[21,37] from already known exact solutions. 

The sine- Gordon equation was the fourth nonlinear partial differential equation whose 
initial-value problem was discovered [2,3] to be solvable by the inverse scattering transform 
method. This method associates (1.1) with the first-order system of ordinary differential 
equations 



where Ux appears in the coefficients as a potential. By exploiting the one-to-one correspon- 
dence between Ux and the corresponding scattering data for (1.4), the inverse scattering 
transform method determines the time evolution u{x, 0) ^ u{x, t) for (1.1) with the help of 
the solutions to the direct and inverse scattering problems for (1.4). The direct scattering 
problem for (1.4) amounts to finding the scattering coefficients (related to the asymptotics 
of scattering solutions to (1.4) as a; — >■ ±oo) when u{x,t) is known for all x. On the other 
hand, the inverse scattering problem consists of finding u{x,t) from an appropriate set of 
scattering data for (1.4). 




(1.4) 
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In this paper we provide several, but equivalent, explicit formulas for exact solutions 
to (1.1). The key idea to obtain such explicit formulas is to express the kernel of a related 
Marchenko integral equation arising in the inverse scattering problem for (1.4) in terms 
of a real triplet {A,B,C) of constant matrices and by using matrix exponentials. Such 
explicit formulas provide a compact and concise way to express our exact solutions, which 
can equivalently be expressed in terms of exponential, trigonometric (sine and cosine), and 
polynomial functions of x and t. This can be done by "unpacking" matrix exponentials in 
our explicit formulas. As the matrix size increases, the unpacked expressions become very 
long. However, such expressions can be evaluated explicitly for any matrix size either by 
hand or by using a symbolic software package such as Mathematica. One of the powerful 
features of our method comes from the fact that our concise and compact explicit solution 
formulas are valid for any matrix size in the matrix exponentials involved. In some other 
available methods, exact solutions are attempted in terms of elementary functions without 
the use of matrix exponentials, and hence exact solutions produced by such other methods 
will be relatively simple and we cannot expect those methods to produce our solutions 
when the matrix size is large. 

Our method is generalizable and applicable to obtain similar explicit formulas for exact 
solutions to other integrable nonlinear partial differential equations, where a Marchenko 
integral equation is used to solve a related inverse scattering problem. We refer the reader 
to [5-7,14,15], where similar ideas are used to obtain explicit formulas for exact solutions to 
the Korteweg-de Vries equation on the half line and to the focusing nonlinear Schrodinger 
equation and its matrix generalizations. 

In our method, with the help of the matrix triplet and matrix exponentials, we easily 
establish the separability of the kernel of the relevant Marchenko integral equation and 
thus solve it exactly by using linear algebra. We then obtain our exact solutions to the 
sine-Gordon equation by a simple integration of the solution to the Marchenko equation. 
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Our method easily handles complications arising from the presence of non-simple poles of 
the transmission coefficient in the related linear system (1.4). Dealing with non-simple 
poles without the use of matrix exponentials is very complicated, and this issue has also 
been a problem [33,44] in solving other integrable nonlinear partial differential equations 
such as the nonlinear Schrodinger equation. 

Our paper is organized as follows. In Section 2 we establish our notation, introduce 
the relevant Marchenko integral equation, and mention how a solution to the sine-Gordon 
equation is obtained from the solution to the Marchenko equation by using the inverse 
scattering transform method. In Section 3 we outline the solution to the Marchenko 
integral equation when its kernel is represented in terms of a triplet of matrices [A, B, C) 
and thus we derive two solution formulas for exact solutions to the sine-Gordon equation. 
In Sections 4 and 5 we show that our explicit solution formulas hold when the input matrix 
triplets come from a larger family; we show that our solution formulas in the more general 
case can be obtained by constructing two auxiliary constant matrices Q and N satisfying 
the respective Lyapunov equations given in Section 4, or equivalently by constructing 
an auxiliary constant matrix P satisfying the Sylvester equation given in Section 5. In 
Section 4 we also show that the matrix triplet {A, B, C) used as input to construct our 
exact solutions to the sine-Gordon equation can be chosen in various equivalent ways and 
we prove that our exact solutions are analytic on the a;t-plane. In Section 5 we also explore 
the relationship between the Lyapunov equations and the Sylvester equation and show how 
their solutions are related to each other in a simple but interesting way. In that section we 
also show that the two solution formulas derived in Section 3 are equivalent. In Section 6 we 
show that those two equivalent solution formulas can be represented in other equivalent 
forms. In Section 7 we evaluate the square of the spatial derivative of our solutions to 
(1.1) by providing some explicit formulas in terms of the matrix triplet (A, S, C), and we 
evaluate the asymptotics of our exact solutions as a; ^ — oo for each fixed t. In Section 8 
we show that the reflection coefficients associated with such solutions are zero, and we 
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also evaluate explicitly the corresponding transmission coefficient. Finally, in Section 9 we 
provide some specific examples of our exact solutions and their snapshots. 

Let us remark on the logarithm and inverse tangent functions we use throughout our 
paper. The log function we use is the principal branch of the complex- valued logarithm 
function and it has its branch cut along the negative real axis while log(l) = 0. The tan~^ 
function we use is the single- valued branch related to the principal branch of the logarithm 
as 

1 1 , fl + iz\ , ^. 1 /i(l-z)\ 

tan-'. = -log (—j, log. = 2Ua.-'(^j, (1.6) 

and its branch cut is {—ioo, —i][J[i, +ioo). For any square matrix M not having eigenvalues 
on that branch cut, we define 

tan-^M:=— / dz [t&n-'^ z]{zl - M)-'^ , (1.6) 

where the contour 7 encircles each eigenvalue of M exactly once in the positive direction 
and avoids the branch cut of tan~^ z. If all eigenvalues of M have modulus less than 1, we 
then have the familiar series expansion 

tan-i M :=M - -M^ + -M^ - -M^ + . . . . 
3 5 7 

For real-valued h{x) that vanishes as a; — >■ +00, the function tein~^{h{x)) always has 
range (— 7r/2,7r/2) when x values are restricted to (xo,-l-oo) for some large xq value; our 
tan~^(/i(a;)) is the continuous extension of that piece from x e {xq, +00) to a; e (—00, +00). 

2. PRELIMINARIES 

In this section we briefly review the scattering and inverse scattering theory for (1.4) 
by introducing the scattering coefficients and a Marchenko integral equation associated 
with (1.4). We assume that u is real valued and that is integrable in x for each fixed t. 
We also mention how a solution to the sine-Gordon equation is obtained from the solution 
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to the Marchenko equation. We refer the reader to the generic references such as [1,4,25,32] 
for the details. 

Two hnearly independent solutions to (1.4) known as the Jost solutions from the left 
and from the right, denoted by ip^X^x^t) and 0(A,a;,t), respectively, are those solutions 
satisfying the respective spatial asymptotics 



V'(A,a;,t) 



0(A,a;,t) = 





—iXx 





+ o(l), X ^ +00, 



+ o(l), X —CO. 



(2.1) 



The scattering coefficients for (1.4), i.e. the transmission coefficient T, the right reflec- 
tion coeflficient R, and the left reflection coefficient L, can be defined through the spatial 
asymptotics 



•i/'(A, X, t) ~ 



L{\t)e- 



-iXa 



4>{x,x,t) = 



T(A) 

^iXx 

T(A) 

g— iAa; 

R{X, t) e'^"' 



+ o(l), X — >■ — oo. 



(2.2) 



T(A) 



+ o(l), X +00, 



where T does not depend on t, and R and L depend on t as 

R{X, t) = R{X, 0) e-^*/(2A) ^ ^) ^ Q) giV(2A) _ 

We recall that a bound state corresponds to a square- integrable solution to (1.4) and such 
solutions can only occur at the poles of the meromorphic extension of T to the upper half 
complex A- plane denoted by C"*". Because u{x, t) is real valued, such poles can occur either 
on the positive imaginary axis, or for each pole not on the positive imaginary axis there 
corresponds a pole symmetrically located with respect to the imaginary axis. Furthermore, 
such poles are not necessarily simple. If Ux is integrable in x for each fixed t and if the 
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transmission coefficient T is continuous for real values of A, it can be proved by elementary- 
means that the number of such poles and their multiplicities are finite. 

With the convention u{x,t) — >■ as a; — >■ +00, it is known that u{x,t) in (1.4) can be 
determined as 

POO 

u{x,t) = -4: drK{r,r,t), (2.3) 

J X 

or equivalently we have 

Ux{x, t) = 4:K{x, X, t), 
where K{x,y,t) is the solution to the Marchenko integral equation 

("OO r-OO 

K{x,y,t) - n{x + y,ty + dv dr K{x,v,t) ^{v + r,t) ^{r + y,t)* = 0, y > x, 

J X J X 

(2.4) 

where the asterisk is used to denote complex conjugation (without taking the matrix 
transpose) and 

1 

n{y, t) = — dX R{X, t) e'^y + J2 e^^^2/-iV(2A,)^ (2.5) 
27r 

provided the poles Xj of the transmission coefficient are all simple. 

The inverse scattering transform procedure can be summarized via the following diagram: 

/ r\\ / r\\ direct scattering at t=0 rr-i/\/-v\r\ 
U{x,0) > Ua;{x,0) > {R{X,0),{Xj,Cj}} 



sine-Gordon solution 



time evolution 



u{x,t) i Ux{x,t) i {i?(A,t),{Aj,Cje-^*/(2^^)}} 

inverse scattering at t 

We note that in general the summation term in (2.5) is much more complicated, 
and the expression we have provided for it in (2.5) is valid only when the transmission 
coefficient T has simple poles at A-, with j = 1, . . . , n on C""". In case of bound states with 
nonsimple poles, it is unknown to us if the norming constants with the appropriate time 
dependence have ever been presented in the literature. Extending our previous results for 



the nonlinear Schrodinger equation [6,7,12,14] to the sine-Gordon equation, it is possible 
to obtain the norming constants with appropriate dependence on the parameter t in the 
most general case, whether the bound-state poles occur on the positive imaginary axis 
or occur pairwise located symmetrically with respect to the positive imaginary axis, and 
whether any such poles are simple or have multiplicities. In fact, in Section 8 we present 
the norming constants and their proper time dependence on t as well as the most general 
form of the summation term that should appear in (2.5). 

When u is real valued, it is known that for real A we have 

R{-X, t) = R{X, t)*, L(-A, t) = L(A, t)*, T(-A) = T(A)*. 

Because u is real valued, as we verify in Section 3, both the kernel 0,{y, t) and the solution 
K[x, y, t) in (2.4) are also real valued, i.e. 

n{y,tr = n{y,t), (2.6) 
Kix,y,tr^K{x,y,t). (2.7) 

3. EXPLICIT SOLUTIONS TO THE SINE-GORDON EQUATION 

Our goal in this section is to obtain some exact solutions to the sine-Gordon equation 
in terms of a triplet of constant matrices. Following the main idea of [6,7] we will replace 
the summation term in (2.5) by a compact expression in terms of a matrix triplet {A, B, C), 
i.e. we will replace 0,{y, t) when i? = by 

n{y,t) = Ce-^y-^~''/^B, (3.1) 

where A, B, C are real and constant matrices of sizes p x p, p x 1, and 1 x p, respectively, 
for some positive integer p. 
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Recall that any rational function /(A) that vanishes as A — > oo in the complex A- plane 
has a matrix realization in terms of three constant matrices A, B, C as 

/(A) = -iC{XI - iA)-^B, (3.2) 

where / is the p x p identity matrix, A has size p x p, B has size p x 1, and C has size 
1 X p for some p. We will refer to {A, B, C) as a matrix triplet of size p. It is possible to 
pad A, B, C with zeros or it may be possible to change them and increase or decrease the 
value of p without changing /(A). The smallest positive integer p yielding /(A) gives us a 
"minimal" realization for /(A), and it is known [8] that a minimal realization is unique up 
to a similarity transformation. Thus, without any loss of generality we can always assume 
that our triplet (A, S, C) corresponds to a minimal realization, and we will refer to such a 
triplet as a minimal triplet. Note that the poles of /(A) correspond to the eigenvalues of 
the matrix {lA). By taking the Fourier transform of both sides of (3.2), where the Fourier 
transform is defined as 

we obtain 

f{y) = Ce-'^yB. (3.3) 

We note that under the similarity transformation (A, C) {S~^AS,S~^B,CS) for 
some invertible matrix S, the quantities /(A) and f{y) remain unchanged. 

Comparing (3.1) and (3.3) we see that they are closely related to each other. As 
mentioned earlier, without loss of any generality we assume that the real triplet {A, B, C) 
in (3.1) corresponds to a minimal realization in (3.2). For the time being, we will also 
assume that all eigenvalues of A in (3.1) have positive real parts. However, in later sections 
we will relax the latter assumption and choose our triplet in a less restrictive way, i.e. in 
the admissible class A defined in Section 4. 

Let us use a dagger to denote the matrix adjoint (complex conjugation and ma- 
trix transpose). Although the adjoint and the transpose are equal to each other for real 
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matrices, we will continue to use the dagger notation even for the matrix transpose of 
real matrices so that we can utilize the previous related results in [5,6] obtained for the 
Zakharov-Shabat system and the nonlinear Schrodinger equation. Since Q appearing in 
(3.1) is a scalar we have 0,^ = fl*; thus, we get 

We note that when Q is given by (3.1), the Marchenko equation is exactly solvable by 
using linear algebra. This follows from the separability property of the kernel, i.e. 

Q{x + y,t) = Ce-^^e-^^-^"'*/25^ (3.4) 

indicating the separability in x and y; thus, (3.4) allows us to try a solution to (2.4) in the 
form 

K{x, y, t) = H{x, t) e-^'y-(^'r't/2c\ (3.5) 
Using (3.5) in (2.4) we get 

H{x, t) r{x, t) = S+e"^^^, (3.6) 

or equivalently 

H{x, t) = 5te-^'^r(a;, t)-\ (3.7) 

where we have defined 

r(a;, t):=I + e-^'--^^'r't/2Q^-2A.-A-^t/2js^^-A^x^ (3_g) 

with the constant p x p matrices Q and N defined as 

POO poo 

Q := dse-'^^'C^Ce-^', N := e-^'^SStg-^^^ (3.9) 

Jo Jo 

It is seen from (3.9) that Q and N are selfadjoint, i.e. 

Q = Q\ N = Nl (3.10) 
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In fact, since the triplet {A, B, C) is real, the matrices Q and N are also real and hence 
they are symmetric matrices. Using (3.7) in (3.5) we obtain 

K{x, y, t) = Ste-^'"r(a;, t)"^ e-A^y-iA^)-'t/2fji ^ ^3_^^) 

or equivalently 

K{x, y, t) = B^F{x, t)-ie-^^(^-^)Ct, (3.12) 

where we have defined 

F{x,t) :^ e^^ +Qe-^N, (3.13) 

with the quantity P defined as 

P := 2Ax + ^ A'H. (3.14) 

From (2.3) and (3.12) we see that 

/•oo 

u{x,t) = -A drB^F{r,t)-^Cl (3.15) 

J X 

The procedure described in (3.4)-(3.15) is exactly the same procedure used in [5,6] 
with the only difi'erence of using A~'^t/2 in the matrix exponential in (3.15) instead of AiA^t 
used in [5,6]. However, such a difference does not affect the solution to the Marchenko 
integral equation at all thanks to the fact that A and A~^ commute with each other. In 
fact, the solution to the Marchenko equation is obtained the same way if one replaces 
A~^ /2 by any function of the matrix A because such a matrix function commutes with A. 

We will later prove that F{x, t) given in (3.13) is invertible on the entire xt-plane and 
that F{x,t)~^ exponentially as a; — > ±oo and hence u{x,t) given in (3.15) is well 
defined on the entire xt- plane. We note that, as a result of (2.6), the solution K{x, y, t) to 
the Marchenko equation (2.4) is real and hence (2.7) is satisfied. Hence, from (2.3) we see 
that u{x,t) is real valued, and by taking the adjoint of both sides of (3.15) we get 

POO 

u{x,t) = -A drC[F{r,t)^]-'^B. (3.16) 

J X 
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Instead of using (2.6) at the last stage, let us instead use it from the very beginning 
when we solve the Marchenko equation (2.4). Replacing Q* by O in the two occurrences 
in (2.4), we can solve (2.4) in a similar way as in (3.4)-(3.15) and obtain 

K{x, y, t) = CE{x, t)-^e-'^^y-''^B, (3.17) 

where we have defined 

E{x,t) := + Pe-'^P, (3.18) 
with P as in (3.14) and the constant matrix P given by 

P:= / dse-^'BCe-^'. (3.19) 
Jo 

Thus, from (2.3) and (3.17) we obtain 

/■oo 

u{x,t) = -Aj drCE{r,t)-^B. (3.20) 

J X 

We will show in Section 5 that the two explicit solutions to the sine-Gordon equation given 
by (3.15) and (3.20) are identical by proving that 

E{x,t) = F{x,t)l (3.21) 



4. EXACT SOLUTIONS USING THE LYAPUNOV EQUATIONS 

In Section 3 we have derived (3.15) and (3.20) by assuming that we start with a real 
minimal triplet (A, B, C) where the eigenvalues of A have positive real parts. In this section 
we show that the explicit formula (3.15) for exact solutions to the sine-Gordon equation 
remains valid if the matrix triplet {A,B,C) used to construct such solutions is chosen in 
a larger class. Starting with a more arbitrary triplet we will construct the matrix F given 
in (3.13), where the auxiliary matrices Q and N are no longer given by (3.9) but obtained 
by uniquely solving the respective Lyapunov equations 

A^Q + QA = C^C, (4.1) 
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AN + NA"^ = BBl 



(4.2) 



Many of the proofs in this section are similar to those obtained earher for the nonhnear 
Schrodinger equation [5,6] and hence we will refer the reader to those references for the 
details of some of the proofs. 

Definition 4.1 We say that the triplet {A,B,C) of size p belongs to the admissible class 
A if the following conditions are met: 

(i) The matrices A, B, and C are all real valued. 

(ii) The triplet (A, B, C) corresponds to a minimal realization for /(A) when that triplet 
is used on the right hand side of (3.2). 

(iii) None of the eigenvalues of A are purely imaginary and no two eigenvalues of A can 
occur symmetrically with respect to the imaginary axis in the complex X-plane. 

We note that, since A is real valued, the condition stated in (iii) is equivalent to the 
condition that zero is not an eigenvalue of A and that no two eigenvalues of A are located 
symmetrically with respect to the origin in the complex plane. Equivalently, (iii) can be 
stated as A and {—A) not having any common eigenvalues. We will say that a triplet is 
admissible if it belongs to the admissible class A. 

Starting with a triplet {A,B,C) in the admissible class A, we will obtain exact solu- 
tions to the sine- Gordon equation as follows: 

(a) Using A, B, C as input, construct the airxiliary matrices Q and N by solving the re- 
spective Lyapunov equations (4.1) and (4.2). As the next theorem shows, the solutions 
to (4.1) and (4.2) are unique and can be obtained as 




(4.3) 



N 




7 



(4.4) 
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where 7 is any positively oriented simple closed contour enclosing all eigenvalues of 
{iA) and leaving out all eigenvalues of {—iA^). If all eigenvalues of A have positive 
real parts, then Q and N can also be evaluated as in (3.9). 

(b) Using the auxiliary matrices Q and N and the triplet (A, 5,(7), form the matrix 
F{x, t) as in (3.13) and obtain the scalar u{x,t) as in (3.15), which becomes a solution 
to (1.1). 

Theorem 4.2 Consider any triplet {A, B, C) belonging to the admissible class A described 
in Definition 4-i- Then: 

(i) The Lyapunov equations and (4.2) are uniquely solvable, and their solutions are 
given by (4-3) and (4-4)} respectively. 

(ii) The constant matrices Q and N given in (4-3) and (4-4)} respectively, are selfadjoint; 
i.e. = Q and ATt = N. In fact, since the triplet (A, B, C) is real, the matrices Q 
and N are also real. Furthermore, both Q and N are invertible. 

(iii) The resulting matrix F{x,t) formed as in (3.13) is real valued and invertible on the 
entire xt-plane, and the function u{x,t) defined in (3.15) is a solution to the sine- 
Gordon equation everywhere on the xt-plane. Moreover, u{x,t) is analytic on the 
entire xt-plane and u^ixji) decays to zero exponentially as x ^ ±00 at each fixed 
t e R. 

PROOF: The proof of (i) follows from Theorem 4.1 of Section 4.1 of [20]. It is directly seen 
from (4.1) that is also a solution whenever Q is a solution, and hence the uniqueness of 
the solution assures Q — . Similarly, as a result of the realness of the triplet {A, B, C), 
one can show that Q* is also a solution to (4.1) and hence Q = Q*. The selfadjointness and 
realness of N are established the same way. The invertibility of Q and N is a result of the 
minimality of the triplet (^4, B, C) and a proof can be found in the proofs of Theorems 3.2 
and 3.3 of [5] by replacing (2.2) of [5] with (3.13) in the current paper, completing the 
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proof of (ii). From (3.13) and (3.14) it is seen that the realness of the triplet {A,B,C) 
and of Q and N imphes the realness of F. The proof of the invertibility of F is similar to 
the proof of Proposition 4.1 (a) of [5] and the rest of the proof of (iii) is obtained as in 
Theorem 3.2 (d) and (e) of [5]. | 

We will say that two triplets {A, B, C) and {A, B, C) are equivalent if they lead to the 
same u{x,t) given in (3.15). The next result shows that two admissible triplets are closely 
related to each other and can always be transformed into each other. 

Theorem 4.3 For any admissible triplet {A,B,C), there corresponds an equivalent ad- 
missible triplet {A, B, C) in such a way that all eigenvalues of A have positive real parts. 

PROOF: The proof is similar to the proof of Theorem 3.2 of [5], where the triplet (Jl, B, C) 
is expressed explicitly when one starts with the triplet {A,B,C). Below we provide the 
explicit formulas of constructing (yl, S, C) by starting with (A, 5,(7); i.e., by providing 
the inverse transformation formulas for those given in [5] . Without loss of any generality, 
we can assume that {A, B, C) has the form 

\bA . . . 

C=[Ci C2], 





"ii 


" 




'Bi' 


A = 






B = 









A2 




B2 



where all eigenvalues of Jii have positive real parts and all eigenvalues of A2 have negative 

real parts, and for some ^ < q < p, the sizes of the matrices Ai, A2, Bi, B2, Ci, C2 are 

q X q, ip — q) X ip — q), q X 1, {p — q) x 1, 1 x q, and 1 x {p — q), respectively. We first 

construct the matrices Q and N by solving the respective Lyapunov equations 

QA + itg = c^c, 



AN + NA^ = BBl 
Writing Q and N in block matrix forms of appropriate sizes as 

Q 

and, for appropriate block matrix sizes, by letting 



'Qi 


Q2 




"iVi 


N2 






iV = 








Q4. 






7V4_ 







" 




'B{ 


A = 




B^ 









A2 




B2 



C2 



(4.5) 



(4.6) 
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Ci = Ci — C2Q4 ^Q^i C2 = ^2(^4 ^ 1 



we obtain 

^i=ii, A2 = -Al Bi = B^-N2N^^B2, B2 = N^'B2, (4.7) 

(4.8) 

yielding {A, B,C) by starting with {A,B,C). | 

When the triplet {A,B,C) is decomposed as in (4.6), let us decompose the corre- 
sponding solutions Q and N to the respective Lyapunov equations (4.1) and (4.2), in an 
analogous manner to (4.5), as 



Q 



'Qi 


Q2 






N2' 






N = 






Qz 








^4. 



(4.9) 



The relationship between (4.5) and (4.9) is summarized in the following theorem. 

Theorem 4.4 Under the transformation (A, B, C) i->- {A, B, C) specified in Theorem 4-3, 
the quantities Q, A^, F, E appearing in (4-.1), (4--2), (3.13), (3.18), respectively, are trans- 
formed as 

{Q,N,F,E)^{Q,N,F,E), 



where 



Q 



Qi — (52Q4 ^Qs — Q2Q4 ^ 

—Q4 ^Qs —Qa ^ 



N 



-n: 



I -Q2Q4 

-Qj' 



/ 
I 

PROOF: The proof can be obtained in a similar manner to the proof of Theorem 3.2 of 

[5] by using 



E = 



I -N2N^ 
-N^^ 



E 



(4.10) 
(4.11) 





'Ai 




A = 


. -4. 





I -N2N^ 
-N^^ 



B, 



C^C 



I 

—Q4 ^Qs —Qa ^ 
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corresponding to the transformation specified in (4.7) and (4.8). | 

As the foUowing theorem shows, for an admissible triplet {A,B,C), there is no loss 
of generality in assuming that all eigenvalues of A have positive real parts and B has a 
special form consisting of zeros and ones. 

Theorem 4.5 For any admissible triplet {A,B,C), there correspond a special admissible 
triplet (A, B, (7), where A is in a Jordan canonical form with each Jordan block containing 
a distinct eigenvalue having a positive real part, the entries of B consist of zeros and ones, 
and C has constant real entries. More specifically, for some appropriate positive integer m 
we have 





\A^ 


• 


■ 






A = 





A2 ■ 





B = 






. 


• 


-^m - 




-Bm- 



(4.12) 



where in the case of a real (positive) eigenvalue Uj of Aj the corresponding blocks are given 
by 



00 -i 



-1 • 



ujj 




-1 




L 













• • Cj2 Cji J , 

1 







ujj J 



Bj 



(4.13) 



(4.14) 



with Aj having size rij x nj, Bj size rij xl, Cj size 1 x n^, and the constant Cjnj is nonzero. 
In the case of complex eigenvalues, which must appear in pairs as aj ± iPj with aj > 0, 
the corresponding blocks are given by 



[ijn, ej 

Aj -I2 . 

Aj . 

A, . 








Uj ■ ■ ■ 




.. 


- 


.. 





.. 





.. A, 


-I2 


.. 


Aj J 



B, 



(4.15) 



(4.16) 
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where and ejs for s = 1, . . . ,nj are real constants with {I'jn^ +^jnj) > 0) -^2 denotes the 
2x2 unit matrix, each column vector Bj has 2nj components, each Aj has size 2nj x 2nj , 
and each 2x2 matrix Aj is defined as 



Aj:= 



(4.17) 

PROOF: The real triplet {A,B,C) can be chosen as described in Section 3 of [7]. | 

5. EXACT SOLUTIONS USING THE SYLVESTER EQUATION 

In Section 3, starting from a minimal triplet (A, S, C) with all eigenvalues of A having 
positive real parts, we have obtained the exact solution formula (3.20) to the sine-Gordon 
equation by constructing the matrix E{x,t) in (3.18) with the help of the auxiliary matrix 
P in (3.19). In this section we show that the explicit formula (3.20) for exact solutions to 
the sine-Gordon equation remains valid if the matrix triplet {A, B, C) used to construct 
such solutions comes from a larger class, namely from the admissible class A specified in 
Definition 4.1. 

Starting with any triplet {A^ B, C) in the admissible class A, we obtain exact solutions 
to the sine- Gordon equation as follows: 

(a) Using A^ B, C as input, construct the auxiliary matrix P by solving the Sylvester 
equation 

AP + PA^ BC. (5.1) 
The unique solution to (5.1) can be obtained as 

P=— [ dX{XI- iA)-^BC{XI + iA)-\ (5.2) 

where 7 is any positively oriented simple closed contour enclosing all eigenvalues of 
(iA) and leaving out all eigenvalues of {—iA). If all eigenvalues of A have positive real 
parts, then P can be evaluated as in (3.19). 
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(b) Using the auxiliary matrix P and the triplet {A^ B, C), form the matrix E{x,t) as in 
(3.18) and then form the scalar u{x, t) as in (3.20). 

Theorem 5.1 Consider any triplet (A, B, C) belonging to the admissible class A described 
in Definition 4-i- Then, the Sylvester equation (5.1) is uniquely solvable, and its solution 
is given by (5.2). Furthermore, that solution is real valued. 

PROOF: The unique solvability of (5.1) is already known [20]. For the benefit of the reader 
we outline the steps below. From (5.1) we get 

-{XI - iA)P + P{\I + lA) = iBC, 

or equivalently 

-P{XI + lA)-^ + {XI - iA)-^P = i{XI - iA)-^BC{XI + lA)-^. (5.3) 

Dividing both sides of (5.3) by {2%) and then integrating along 7, and using 

J- [ dX {XI - iA)-^ = /, ^ I dX {XI + lA)-^ = 0, 

we obtain (5.2) as the unique solution to (5.1). Since the admissible triplet {A,B,C) is 
real, by taking complex conjugate of both sides of (5.1) we see that P* also solves (5.1). 
From the uniqueness of the solution to (5.1), it then follows the P* = P. | 

Next we show that, for any triplet {A, B, C) in our admissible class A, there is a close 
relationship between the matrix P given in (5.2) and the matrices Q and N appearing in 
(4.3) and (4.4), respectively. 

Theorem 5.2 Let the triplet {A,B,C) of size p belong to the admissible class specified in 
Definition Then the solution P to the Sylvester equation (5.1) and the solutions Q 
and N to the respective Lyapunov equations (4-1) o-nd (4-2) satisfy 

NQ = P^. (5.4) 
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PROOF: Note that (5.4) is valid when the matrix A in the triplet is diagonal. To see this, 
note that the use of the triplet (A, S, C) with 

■bi 



A = diag{ai, • • • ,ap}, 

in (4.1), (4.2), and (5.1) yields 

bjCk 



B 



C=[ci 



-P J 5 



aj + ak 



Qjk — 



-bp. 



CjCk 

aj + ak 



Njk = 



bjbk 
aj + ak 



where the subscript jk denotes the (j, k) entry of the relevant matrix. Hence, 



iNQ)jk = j2 



bjbsCgCk 



bjCsbsCk 



{aj + as){as + ak)' ^"^^ ^ (oj + Os)(os + Ofe) ' 

establishing (5.4). Next, let us assume that A is not diagonal but diagonalizable through 
a real-valued invertible matrix S so that A — S~^AS and A is diagonal. Then, under the 
transformation 

{A, B, C) ^ {A, B, C) = {S-'^AS, S'^B, CS), 



we get 



{Q,N,P) ^ {Q,N, P) = {S^QS, S-^N{S^)-\S-^PS), 



where Q, N, and P satisfy (4.1), (4.2), and (5.1), respectively, when {A,B,C) is replaced 
with (A, C) in those three equations. We note that (A, B, C) is an admissible triplet 
when {A,B,C) is admissible because the eigenvalues of A and A coincide. Since A is 
diagonal, we already have NQ = P^, which easily reduces to NQ = P^ given in (5.4). In 
case A is not diagonalizable, we proceed as follows. There exists a sequence of admissible 
triplets {Ak, B, C) converging to {A, B, C) as k ^ +oo such that each Ak is diagonalizable. 
Let the triplet {Qki ^ki Pk) correspond to the solutions to (4.1), (4.2), and (5.1), respec- 
tively, when {A, B, C) is replaced with {Ak, B, C) in those three equations. We then have 
NkQk = -Pfc 5 and hence {Qk, Nk, Pk) {Q, N, P) yields (5.4). Note that we have used the 
stability of solutions to (4.1), (4.2), and (5.1). In fact, that stability directly follows from 
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the unique solvability of the matrix equations (4.1), (4.2), (5.1) and the fact that their 
unique solvability is preserved under a small perturbation of A. | 

Theorem 5.3 Let the triplet (A, B, C) belong to the admissible class specified in Definition 
4-1- Then, the solution P to the Sylvester equation (5.1) and the solutions Q and N to 
the respective Lyapunov equations (4.1) and (4.2) satisfy 

N{A^yQ = PA^P, j = 0,±1,±2,.... (5.5) 

PROOF: Under the transformation 

{A, B, C) ^ {A, B, C) = {A, A^B, C), 

we get 

(g, N, P) ^ {Q, N, P) = (Q, A^N{A^y, A^P), 

where Q, N, and P satisfy (4.1), (4.2), and (5.1), respectively, when {A,B,C) is replaced 
with [A, B, C) in those three equations. Since {A, B, C) is also admissible, (5.4) implies 
that NQ = P^, which yields (5.5) after a minor simplification. | 

Next, given any admissible triplet {A, B, C), we prove that the corresponding solution 
P to (5.1) is invertible and that the matrix E{x,t) given in (3.18) is invertible and that 
(3.21) holds everywhere on the a;t-plane. 

Theorem 5.4 Let the triplet (A, B, C) belong to the admissible class specified in Definition 
4.1, and let the matrices Q, N, P be the corresponding solutions to (4-1), (4-2), and (5.2), 
respectively. Then: 

(i) The matrix P is invertible. 

(ii) The matrices F and E given in (3.13) and (3.18), respectively, are real valued and 
satisfy (3.21). 

(iii) The matrix E{x,t) is invertible on the entire xt-plane. 
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PROOF: The invertibility of P follows from (5.4) and the fact that both Q and N are 
invertible, as stated in Theorem 4.2 (ii); thus, (i) is established. To prove (ii) we proceed 
as follows. The real-valuedness of F has already been established in Theorem 4.2 (iii). 
From (3.18) it is seen that the real-valuedness of the triplet {A^B^C) and of P implies 
that E is real valued. From (3.13), (3.14), and (3.18) we see that (3.21) holds if and only 
if we have 

Ne-^^Q = Pe-^P, (5.6) 

where we have already used N''' = N and = Q, as established in Theorem 4.2 (ii). Since 
(5.5) implies 

N{-(3^yQ = P{-(3yP, j = 0,1,2,..., 

we see that (5.6) holds. Having established (3.21), the invertibility of E{x,t) on the entire 
a;t-plane follows from the invertibility of F{x, t), which has been established in Theorem 4.2 
(iii). I 

Next, we show that the explicit formulas (3.15), (3.16), and (3.20) are all equivalent 
to each other. 

Theorem 5.5 Consider any triplet (A, B, C) belonging to the admissible class A described 
in Definition 4-i- Then: 

(i) The explicit formulas (3.15), (3.16), and (3.20) yield equivalent exact solutions to the 
sine-Gordon equation (1-1) everywhere on the entire xt-plane. 

(ii) The equivalent solution u{x,t) given in (3.15), (3.16), and (3.20) is analytic on the 
entire xt-plane, and Ux{x,t) decays to zero exponentially as x ^ ±oo at each fixed 
t e R. 

PROOF: Because u{x, t) is real and scalar valued, we already have the equivalence of (3.15) 
and (3.16). The equivalence of (3.16) and (3.20) follows from (3.21). We then have (ii) as 
a consequence of Theorem 4.2 (iii). | 
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6. FURTHER EQUIVALENT FORMS FOR EXACT SOLUTIONS 

In Theorem 5.5 we have shown that the exact solutions given by the exphcit formulas 
(3.15), (3.16), and (3.20) are equivalent. In this section we show that our exact solutions 
can be written in various other equivalent forms. We first present two propositions that 
will be useful in later sections. 

Proposition 6.1 If {A, B,C) is admissible, then the quantities and E~^, appearing 
in (3.13) and (3.18), respectively, vanish exponentially as x ^ ±oo. 

PROOF: It is sufficient to give the proof when the eigenvalues of A have all positive real 
parts because, as seen from (4.10) and (4.11), the same result also holds when some or all 
eigenvalues of A have negative real parts. When the eigenvalues of A have positive real 
parts, from (3.13) we get 

^-1 = e-^'/2[/ + e-^'/2ge-^iVe-^'/2]-le-^'/^ (6.1) 

where the invertibility of Q and N is guaranteed by Theorem 4.2 (ii). Hence, (6.1) implies 
that — > exponentially as a; — > +oo. From (3.21) and the realness of E and F we 
also get exponentially as a; — >■ +oo. To obtain the asymptotics as a; — )■ — oo, we 

proceed as follows. From (3.13) we obtain 

and hence 

and thus F~^ — > exponentially as a; — > — oo. From (3.21) and the realness of E and F 
we also get E~^ — > exponentially as a; — > — oo. | 

Proposition 6.2 The quantity E{x,t) defined in (3.18) satisfies 

E^ = 2AE - 2BCe-^P, Ee'^P = Pe'^E, e^P'^E = P + e^P'^e^. (6.2) 

24 



If {A, B, C) is admissible and all eigenvalues of A have positive real parts, then E ^Pe ^ — > 
exponentially as x ^ —oo. 

PROOF: We obtain the first equafity (6.2) by taking the a;-derivative of (3.18) and by 
using (5.1). The second equality can be verified directly by using (3.18) in it. The third 
equality is obtained by a direct premultiplication from (3.18). The limit as a; ^ — oo is 
seen from the last equality in (6.2) with the help of (3.14). | 

Let us start with a triplet {A, B, C) of size p belonging to the admissible class specified 
in Definition 4.1. Letting 

M(x,t) :=e-^/2pe-^/^ (6.3) 

where ^ as in (3.14) and P is the unique solution to the Sylvester equation (5.1), we can 
write (3.18) also as 

£;(a;,t) = e^/2Ae^/2, (6.4) 

where we have defined 

K{x,t):=I+[M{x,t)f. (6.5) 
Using (5.1) in (6.3), we see that the x-derivative of M(a;, t) is given by 

M^{x,t) = -e-^''^BCe-^'^. (6.6) 

Proposition 6.3 The eigenvalues of the matrix M defined in (6.3) cannot occur on the 
imaginary axis in the complex plane. Furthermore, the matrices (/ — iM) and (/ + iM) 
are invertible on the entire xt-plane. 

PROOF: From (6.4) and (6.5) we see that 

(/ - iM){I + iM) = e-^/2^e-^/2^ 

and by Theorem 5.4 (iii) the matrix E is invertible on the entire xt-plane. Thus, both 
{I — iM) and (I + iM) are invertible, and consequently M cannot have eigenvalues ±i. For 
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any real, nonzero c, consider the transformation {A, B, C) i-> {A^ cB, cC) of an admissible 
triple {A, B, C). The resulting triple is also admissible, and as seen from (5.1) and (6.3) we 
have (P, M, / + M^) i-^ (c^P, c^M, / + c^M^). Thus, M cannot have any purely imaginary 
eigenvalues. Since P is known to be invertible by Theorem 5.4 (i), as seen from (6.3) the 
matrix M is invertible on the entire a;t-plane and hence cannot have zero as its eigenvalue. | 

Theorem 6.4 The solution to the sine-Gordon equation given in the equivalent forms 
(3.15), (3.16), and (3.20) can also he written as 

w(a;,t) = -4Tr[tan-^M(x,t)], (6.7) 

.(.,t)^2dog( ^^^|^_^^|^;^jj j, (6.8) 

nix t) - 4tan-i (^ det(/ + iM(x, t)) - det(/ - .M(x, t)) \ 

u[x, t) - tan ^^^^^ ^ .^^^^ ^ ^^^^^ _ .^^^^ J , (b.yj 

where M is the matrix defined in (6.3) and Tr denotes the matrix trace (the sum of diagonal 
entries). 

PROOF: Let us note that the equivalence of (6.8) and (6.9) follows from the second equality 
in (1.5) by using z = det(/ + iM)/ det(/ — iM) there. To show the equivalence of (6.7) 
and (6.8), we use the matrix identity 

tan"^ M = -Jt log ((/ + iM){I - iM)'^) , 

which is closely related to the first identity in (1.5), and the matrix identity 

Tr[log2;] = logdetz, 

with the invertible matrix z = {I -\- iM){I — iM)~^. Thus, we have established the equiv- 
alence of (6.7), (6.8), and (6.9). We will complete the proof by showing that (3.20) is 
equivalent to (6.7). Using the fact that for any m x n matrix a and any n x m matrix 7 
we have 

Tr[a7] = Tr[7a], (6.10) 
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from (6.4)-(6.6) we get 

-4CE-^B ^ 4Tr[M^{I + M^)-^]. (6.11) 

By Proposition 6.1 we know that vanishes exponentially as x ^ +00. Hence, with 
the help of (6.11) we see that we can write (3.20) as 



u{x,t) = 4Tr 



■ /"OO 

/ drMr{r,t)[I + M{r,t)M{r,t)]-^ 

J X 



which yields (6.7). | 



Theorem 6.5 The solution to the sine-Gordon equation given in the equivalent forms 
(3.15), (3.16), (3.20), (6.7)-(6.9) can also he written as 

p 

u{x,t) = -4^ tan-'^ Kj{x,t), (6.12) 

where the scalar functions Kj{x,t) correspond to the eigenvalues of the matrix M{x,t) 
defined in (6.3) and the repeated eigenvalues are allowed in the summation. 

PROOF: At a fixed {x, t)-value, using the matrix identity 

p 

TT[Mix, ty] = (x, t)Y, 5 = 1,2,3,..., 

for large \z\ values in the complex 2;-plane we obtain 

00 00 p p 

s=0 s=0 j=l j=l 

where we dropped the arguments of M and kj for simplicity. Choosing the contour 7 as 
in (1.6) so that each eigenvalue Kj{x,t) is encircled exactly once in the positive direction, 
we can extend (6.13) to 2; e 7 by an analytic continuation with respect to z. Using (6.12) 
in (1.6), we then obtain 

j [tan"^ z] Ti[{zl - = ^^^^ J dz [tan"^ z] [z - Kj)~^, 
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or equivalently 

p 

Tr[tan"^ M(a;, t)] = tan"^ ^^(a;, t), 

which yields (6.12) in view of (6.7). | 

Let us note that the equivalence of (6.7)-(6.9), and (6.12) implies that one can replace 
M by its Jordan canonical form in any of those four expressions without changing the 
value of u{x,t). This follows from the fact that u{x,t) in (6.8) remains unchanged if M 
is replaced by its Jordan canonical form and is confirmed in (6.12) by the fact that the 
eigenvalues remain unchanged under a similarity transformation on a matrix. 

The next result shows that we can write our explicit solution given in (6.12) yet 
another equivalent form, which is expressed in terms of the coefficients in the characteristic 
polynomial of the matrix M[x, t) given in (6.8). Let that characteristic polynomial be given 

by 

p p 

det {zl - M{x,t)) = n [-^ ~ i^ji^^t)] = J](-l)Vj(a;,t) z^'^ 

where the coefficients aj (x, t) can be written in terms of the eigenvalues Kj (x, t) as 

p p 

(70 = 1, C7i=^Kj, (72= ^ KjKk, ... , CTp = Ki ■ ■ ■ Kp, (6.14) 

where we have dropped the arguments and have written kj and aj for Kj{x, t) and crj{x, t), 
respectively, for simplicity. 

Theorem 6.6 The solution to the sine-Gordon equation given in the equivalent forms 
(3.15), (3.16), (3.20), (6.7)-(6.9), and (6.12) can also he written as 



u{x, t) — —4 tan' 



/L(P-1)/2J 

J2 (-l)V2s+l(x,0 
s=0 



v 



LP/2J 

J](-l)V2.+l(x,t) 

s=0 



(6.15) 



/ 



where \J\ denotes the greatest integer function of j and the quantities aj{x,t) are those 
given in (6. 14). 
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PROOF: When p = 2, by letting r]j := tan ^ K,j{x,t) and using the addition formula for 
the tangent function, we obtain 

tanr/i +tanr72 ki + K2 (^1 

tan(?7i + r]2) = 7- rr- r = = , (6.16) 

1 - (tanr7i)(tan?72) 1 - (Tq - a2 

and hence the application of the inverse tangent function on both sides of (6.16) yields 
(6.15). For larger values of p, we proceed by induction with respect to p and by the further 
use of the addition formula for the tangent function. | 

7. FURTHER PROPERTIES OF OUR EXACT SOLUTIONS 

In this section we derive an explicit expression, in terms of a matrix triplet, for the 
square of the spatial derivative of our exact solutions to (1.1) and analyze further properties 
of such solutions. 

Theorem 7.1 If {A,B,C) is admissible, then the solution to the sine-Gordon equation 
given in the equivalent forms (3.15), (3.16), (3.20), (6.7)-(6.9), (6.12), and (6.15) satisfy 

[u,{x,t)f = Tr[(A-iA,),] = Tt[{E-^E,),] = Tr[(F-iF,),], (7.1) 

where A, and F are the quantities appearing in (6.5), (3.18), and (3.13), respectively. 
Consequently, we have 

^'log(detA(x,t)) _ dHog{dei E{x,t)) _ dHog{det F{x,t)) 

PROOF: Let us use the notation of Theorems 4.3 and 4.4 and use a tilde to denote the 
quantities associated with the triplet (A, i?,C), where some or all eigenvalues of A have 
negative real parts. Because of the equivalence stated in Theorems 4.3 and 4.4, we can 
convert the starting triplet (A, S, C) into an admissible triplet (A, S, C) where the matrix 
A has eigenvalues with positive real parts. We will first establish (7.1) and (7.2) for the 
quantities associated with the triplet (A, S, C) and then show that those formulas remain 
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valid when we use {A, B, C) as the input triplet. We exploit the connection between (1.4) 
and the Zakharov-Shabat system given in (2.1) of [6], where q = —iUx/2 and u is real 
valued. From (2.4) and (2.10) of [6] we see that 

Mx,t)f = 8^^i|^, dr[ur{r,t)r = -8G{x,x,t), (7.3) 

where we have 

poo 

G{x, y,t) = - drK{x, r, t)* n{r + y, t)*, (7.4) 

J X 

with K[x, y, t) given in the equivalent forms (3.12) or (3.17), and Q{r + y, t) given in (3.4). 
Since our triplet {A, B, C) is real, both K and Q, are real valued and we can ignore the 
complex conjugations in the integrand in (7.4). Thus, we get 

/•OO 

G{x,y,t) = -CE{x,t)-^ cZre-^(''-^)SCe-^(^+^)-^"'*/^S, (7.5) 

J X 

which is evaluated with the help of (3.19) as 

G(x, y, t) = -CE{x, t)-^Pe-f^e-^^y-^^B, (7.6) 

where /3 is the quantity in (3.14). Omitting the arguments {x,t) and using (7.6) in (7.3) 
we get 

ul = -8[CE-^Pe-^B]x. (7.7) 
Using (6.3) and (6.4) in (7.7) we obtain 

ul = -8[Ce-^/2A-iMe-^/25]^, (7.8) 

where M is the quantity defined in (6.3). With the help of (6.10) we write (7.8) as 

ul = -8TT[e-^/^BCe-^/''A-^M],, 

or equivalently, after using (6.6), we get 

ul = -8IV[M^A-iM]^, 
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Using (6.5) and the fact that M and A ^ commute, we obtain the first equahty in (7.1). 
With the help of (6.4) we obtain 

= -AA -AA + e-^/^E^e-'^/^ 

A-^A^ = -A~^AA -A + A~^e~^/^E^e-^/^, (7.9) 
and hence using (6.4) and (6.10), from (7.9) we obtain 

IV[A-iA^] = -2Tv[A] + Tv[E-^E^], (7.10) 

estabhshing the second equahty in (7.1). With the help of (3.21) and the fact that E and 
F are real valued, we establish the third equality in (7.1). Using the matrix identity 

^ r -1 1 1 5detQ! 51og(detQ:) 
det a ox ox 

we write (7.1) in the equivalent form of (7.2). Now, if we use {A, B, C) instead of (A, S, C), 
we see from (4.11) that, for some constant invertible matrices Y and Z, we have 

E = YEZ, E-^ = Z-^E-^y-^, Ea; = YE^Z, (7.11) 

and hence, with the help of (6.10) and (7.11) we get 

Tv[E-^Ea,] = TT[E-^E^]. (7.12) 

Similarly, (4.10) yields 

F = Z^FY\ P-^ = (rt)-ii?-i(zt)-i, = Z^FaY\ 

which yields 

IV[F-iF^] -Tr[F-iF,]. (7.13) 
Note that from (7.10) and (7.12) we get 



Tr[A-^A^] + 2Tr[i] = TrfA-^A^;] + 2Tr[yl]. 
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(7.14) 



Thus, by taking the ^-derivatives of both sides in (7.12), (7.13), and (7.14), we estabhsh 
(7.1) and (7.2) without any restriction on the sign of the real parts of the eigenvalues of 



Next, we show that the proof of Theorem 7.1 can be obtained directly without using 
(7.3)-(7.6). For this purpose, it is sufficient for us to show that (7.7) can directly be derived 
from (3.20). 

Proposition 7.2 The equality in (3.20) implies (7.7), i.e. we have 



PROOF: We directly evaluate the left hand side of (7.15) by taking the x-derivative of 



(6.2), and we obtain the right hand side in (7.15). | 

The next result shows that tt(— oo,t) must be an integer multiple of 27r. In fact, we 
have w(— oo, t) = 27rj, where j G {— —p + 1, . . . , 0, . . . ,p — l,p}, with p denoting the size 
of the triplet (A, S, C) used to construct our exact solutions. 

Theorem 7.3 // (A, S, C) is admissible and the eigenvalues of A have positive real parts, 
then the solution to the sine- Gordon equation given in the equivalent forms (3.15), (3.16), 



A. I 



(7.15) 



E ^Pe ^ . We simplify the resulting expression by using the first two equalities given in 



(3.20), (6.7)-(6.9),(6.12), and (6.15) satisfies 




(7.16) 



and u{x,t) converges to an integer multiple of {2n) as x ^ —oo. 



PROOF: From (7.6) and the second equation in (7.3) we see that 




and hence with the help of Proposition 6.1, (5.1), and (6.10) we get 



oo 



dr[ur{r,t)f = 8CP-^B = 8Tr[BCP-^] = 8Tr[{AP + PA)p-^] = 16Tr[yl], 



— oo 
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yielding (7.16). By taking the time derivative of both sides of (7.16), we get 

/oo /-oo 
dr Ur{r,t) Urt{r,t) = / drur{r,t) sin(tt(r, t)) = cos(tt(— oo, t)) — cos(tt(+oo, t)), 
-oo J —oo 

which proves that u{—oo,t) is an integer multiple of (27r) because we use the convention 
that ^(+00,*) = 0. I 



8. TRANSMISSION COEFFICIENT AND NORMING CONSTANTS 

In this section we show that our exact solutions given in equivalent forms (3.15), (3.16), 
(3.20), (6.7)-(6.9), (6.12), and (6.15) correspond to zero reflection coefficients in (1.4), 
we evaluate that corresponding Jost solution explicitly in terms of our triplet {A,B,C), 
determine the transmission coefficient explicitly in terms of the matrix A, and we also 
relate our triplet to the norming constants for (1.4) and to their time evolutions. As we 
have seen in Section 4 there is no loss of generality in choosing our triplet in the special 
form specified in Theorem 4.5, and hence in this section we will assume that (A, S, C) has 
the particular form given in (4.12)-(4.17). 

The Jost solution i(j{X,x,t) satisfying the asymptotics (2.1) is given, as in (2.9) of [6], 

by 






/•oo 


K{x,y,t) 




+ dy 


J X 


_G{x,y,t) _ 



(8.1) 



■^(A, x, t) = e 



iXx 



(8.2) 



where K{x,y,t) and G{x,y,t) are the quantities in (3.17) and (7.6), respectively. Using 
(3.17) and (7.6) in (8.1) we obtain 

iCE{x,t)-\XI + iA)-^B 
l-iCE{x,t)-^Pe->^{XI + iA)-^B _ 

where E and /3 are the quantities appearing in (3.18) and (3.14), respectively. With the 
help of Propositions 5.1 and 5.2, by taking the limit of (8.2) as x — ?> — oo and by comparing 
the result with (2.2), we see that L{X,t) = 0, and hence [1,4,32] also R{X,t) = 0, and 

1 



T(A) 



= l-iCP-\XI + iA)-^B. 



(8.3) 
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Using (5.1), with the help of Proposition 4.2 of [6], we can invert (8.3) to get 

T(A) = 1 + iC{XI - iA)-^p-^B. 



(8.4) 



By using Proposition 4.3 of [6] and (5.1), we can write (8.4) as the ratio of two determinants 



as 



_ det(A/ + zA) 
^^^^ - det(A/-a)- ^^-^^ 

Having determined the transmission coefficient T in terms of the matrix A appearing 
in (4.12), let us clarify the relationship between A and the poles and zeros of T in C"*". 
From (8.5) we see that the zeros and poles of T occur exactly at the eigenvalues of {—iA) 
and of (lA), respectively, and that the poles of T occur either on the positive imaginary 
axis or they are pairwise symmetrically located with respect to the imaginary axis in C"*". 
A comparison of T given in (8.5) with Aj given in (4.14) shows that a bound-state pole 
Xj of T located on the positive imaginary axis is related to the eigenvalue Uj of Aj in the 
form Xj = iuij. A comparison of the poles of T given in (8.5) with Aj of (4.16) reveals 
the relationship between the poles off the imaginary axis and the real constants aj and Pj 
appearing in Aj] namely, the pair of bound-state poles of T symmetrically located with 
respect to the imaginary axis in C"'' occur at A = Aj and A = ~A*, where 

Aj = —Pj + iaj, —Xj — Pj + iaj. (8-6) 

Having clarified the relationship between the matrix A appearing in (4.12) and the 
bound-state poles in C"*" of the transmission coefficient T, let us now discuss the relationship 
between the bound-state norming constants and the row vector C appearing in (4.12). 
In case of nonsimple bound-state poles of T, the bound state norming constants can be 
introduced [6,12] in such a way that the generalization from the simple to the nonsimple 
bound states is the most natural. The summation term in (2.5) assumes that there are 
n simple bound-state poles of T at A = Aj with the norming constants Cje"**/^^'^^^. Let 
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us now generalize it to the case where each bound-state pole has multiplicity Uj, i.e. 
when there are rij linearly independent solutions to (1.4) for X = Xj. The most natural 
generalization is obtained by the association 

Cj Cj, Xj iAj, 1 Bj, 

where Aj, Bj, Cj are the matrices appearing in (4.12). The summation term (2.5) then 
generalizes to one of the equivalent terms given in the set of equalities 

j^l j = l s=l 

where {A,B,C) is the special triplet appearing in (4.12) and 9js{t) are the norming con- 
stants associated with the eigenvalue with multiplicity rij. 

From (8.7) we observe the relationship between the bound-state norming constants 
Ojs{t) and the vectors Cj appearing in (4.12). If A^ occurs on the positive imaginary axis, 
then we see that ^js(O) is the same as Cjs appearing in (4.13) and hence the time evolution 
Ojs{0) 1-^ ^ja{t) is governed by 

[Oju^t) ■■■ Oj^it)] = [ejn,{0) ■■■ ^,■l(0)]e-^."*/^ (8.8) 

where Aj is the matrix obtained as in (4.14) by using Uj = —iXj there. We note that the 
norming constants Ojsit) are all real (positive, negative, or zero) with the understanding 
that Cjnj (t) 0. 

Because of the real valuedness stated in (2.6), if the bound-state pole A^ of T occurring 
off the positive imaginary axis has Ojs{t) as the norming constants, then the bound-state 
pole occurring at {—Xj) has Ojs{t)* as the norming constants. In this case (8.6) holds, and 
a comparison of (8.8) with (4.15) and (4.16) reveals that the contribution from the pair Xj 
and (—A*) is given by one of the equivalent forms 
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where {Aj,Bj,Cj) is the real triplet of size 2nj appearing in (4.15) and (4.16). Thus, 
we see that the real constants ejs and jjs appearing in (4.15) are related to the real and 
imaginary parts of the norming constants Ojs{t) as 

ejs = M^js (0)] , 7^-. = -Ini[^,-, (0)] . 
Defining the real 1 x {2nj) vector 

we obtain the time evolution 9js{0) ^jsit) as 

where Aj is the {Irij) x (2nj) matrix appearing in (4.16). 

Let us note that, by using (8.8), we can describe the time evolution of the (complex) 
norming constants Ojs{t) for s = 1, . . . , corresponding to the complex \j given in (8.6) 
by simply replacing the real matrix Aj of size Uj x Uj given in (4.14) with a complex- valued 
Aj of the same size. That complex Aj is simply obtained by replacing ujj in (4.14) by the 
complex quantity {—i\j). In that case, the time evolution of the norming constants 9js{t)* 
for s = 1, . . . ,nj corresponding to the complex —A* given in (8.6) is simply obtained by 
taking the complex conjugate of both sides of (8.8). 

In short, in the most general case the summation term in (2.5) is given by the expres- 
sion Ce~^^~^~^*/^, where the triplet {A,B,C) has the form (4.12). 

9. EXAMPLES 

Example 9.1 The triplet {A,B,C) with 

A=[a], 5 = [1], C=[c], 
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where a > and c 7^ 0, through the use of (3.19) and (6.3), yields 



C 

2a 



M 



^-lax-t I {la 

2a 



and hence from (6.7) we get 

u{x,t) = -4tan-i (^^ e-2«-*/(2«)) 



(9.1) 



If c > 0, the solution in (9.1) is known as a "kink" [25]; it moves to the left with speed 
l/(4a^) and u{x,t) — >■ — 27r as a; — >■ —00. If c < 0, the solution in (9.1) is known as an 
"antikink" [25]; it moves to the left with speed l/(4a^) and u{x, t) — >■ 27r as a; — >■ —00. 

Example 9.2 The triplet {A,B,C) with 



A = 



a b 
—b a 



B = 



C =[C2 Ci ] , 

where a > 0, 6 7^ 0, and C2 ^ 0, through the use of (3.19), (6.3), and (6.9) yields 



u{x,t) = —4 tan — — , 
\ den / 



(9.2) 



where 



num := 8a^e"^+ [(aci — 6C2) cos(6C-) — (bci + ac2) sin(6C_)] , 

t 



den := b^{cl + cl) + 160^(0^ + b^)e^''^+ , C± := 2a; ± 



2{a? + b'^)' 

The solution in (9.2) corresponds to a "breather" [25] and u{x^ t) ^ as a; ^ —00. For 
example, the choice a = 1, 6 = 2, ci = 2, C2 = 1 simplifies (9.2) to 

.1 /2e2^+*/i°sin(4a; -t/5)' 



u{x^ t) = 4 tan 
Example 9.3 The triplet (A, S, C) with 



V 1 + 46"^^+*/^ 



ai 


" 




"1" 





a2 


1 



C = [ Cl C2 ] , 
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where ai and 02 are distinct positive constants, and ci and C2 are real nonzero constants, 
by proceeding the same way as in the previous example, yields (9.2) with 

num := 2{a, + as)" (aiCae^^^^+^/^^^i) + a2Cie2"^-+*/(2«2)) , 

den := -(ai - 02)^01 C2 + 4aia2(ai + a2)2e^«i+«2)(2x+t/(2aia2))_ ^g g^ 

If (C1C2) < then the quantity in (9.3) never becomes zero; the corresponding solution is 
known as a "soliton-antisoliton" [25] interaction. On the other hand, if (C1C2) > then the 
quantity in (9.3) becomes zero on a curve on the a;t- plane and the corresponding solution 
is known as a "soliton-soliton" [25] interaction. For example, the choice a = 1, 6 = 2, 
ci = ±1, C2 = ^1 yields 

.1 /18e2-+*/2 _36e4-+*/4)- 



u(x,t) = ±4 tan' 



1 + 72e6^+3t/4 



= 2, ci = ±1, 



with u{x,t) — > as a; — > —00. On the other hand, the choice a = 1, b 
C2 = ±1 yields the solution 

uix t) - ^4 tan- ^18e^!!!{!±36e^\ 
u{x,t) - ^4 tan 1^ _^ ^ 72e6^+3V4 J ' 

with u{x, t) ^Att as X ^ —00. 
Example 9.4 The triplet {A, B, C) with 



C = [ C3 C2 Ci ] , 



where a > 0, and ci, C2, C3 are real constants with C3 7^ 0, by proceeding the same way as 
in the previous example, yields u{x,t) in the form of (9.2), where 

num := c^e"^'*^-*/" + 32^?, den := 4ae-2«^-*/(2«) [128a^e^«*+*/" + hi + 7*2] , 





a 


-1 


" 




"0" 


A = 





a 


-1 


B = 













a 




1 



g := {8a'^ci+8a'^C2+8a^C3)-{4:a^C2+8ac3)t+C3r + {16a^C2 + 16a'^C3)x-8a^C3xt+16a'^C3X^, 
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hi := (Sa^cl - Sa^ciCs + 16a^C2C3 + lAa^cj) - (4a^C2C3 + 4:acl)t, 
h2 := cit^ + (16a^C2C3 + 32a^cl)x - Sa^cltx + IGa^clx^. 

The choice a — 1, ci — —1, C2 = —1, cs = —2 yields 

_i / e-^^-* + 8(16- lot + + 24a; -8fa+16a;^) \ 
tan (^2e-2»'-*/2[32e4»'+* + 20 - 6t + _ gta; + 40a; + 16x^] J ' 

with t) — )■ 27r as X — )■ — cxo. On the other hand, the choice a = 1, ci = 0, C2 = 0, C3 = 1 
yields 

, , , _i / e-^^-* + 32(8 -8t + t'^ + 16a; - 8ta; + 16a;2) \ 

uix.t) = —4 tan — -j^-z r— 7- — , 

^ ' ^ \4e-^^-^/^[128e^^+* + 14 - 4t + _ gtx + 32x + 16a;2] J ' 

with u{x, t) — >■ — 27r as a; — >■ —00. 
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